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The violation of the Fluctuation-Dissipation Theorem (FDT) in a two-dimensional Ising model 
with both ferromagnetic exchange and antiferromagnetic dipolar interactions is established and 
investigated via Monte Carlo simulations. Through the computation of the autocorrelation C(t + 
t w ,t w ) and the integrated response (susceptibility) functions we obtain the FDT violation factor 
X(t + t w ,t w ) for different values of the temperature, the waiting time t w and the quotient 8 = Jo/ Jd, 
Jo and Jd being the strength of exchange and dipolar interactions respectively. For positive values of 
5 this system develops a striped phase at low temperatures, in which the non-equilibrium dynamics 
presents two different regimes. Our results show that such different regimes are not reflected in the 
FDT violation factor, where X goes always to zero for high values of t w in the aging regime, a result 
that appears in domain growth processes in non-frustrated ordered systems. 

PACS numbers: 75.40. Gb, 75.40.Mg, 75.10.Hk 



The competition between long-range antiferromagnetic dipolar interactions and short-range ferromagnetic exchange 
interactions can give rise to a variety of unusual and interesting macroscopic phenomena. Recent works in two 
dimensional uniaxial spin systems, where the spins are oriented perpendicular to the lattice and coupled with these 
kind of interactions, have shown a very rich phenomenological scenario concerning both its equilibrium statistical 
mechanicsalj and non-equilibrium dynamical properties!! Moreover, recent results have shown some similarities 
between the non-equilibrium dynamical properties of these kind of ordered systems and that of glassy systemsu. 

Magnetization processes in these kind of systems are of interest due to aspects related to information storage in 
ultrathin ferromagnetic films. Moreover, there are several contexts in which a short-ranged tendency to order is 
perturbated by a long-range frustrating interaction. Among others, model systems of this type has been proposed to 
study avoided phase transitions in supercooled liquidaJ and charge density waves in doped antiferromagnetsou 

The above mentioned systems can be described by an Ising like Hamiltonian of the type 

H=-6j2a i a j+ j:^i (1) 

<ij> « 

where the spin variable Ui = ±1 is located at the site i of a square lattice, the sum 53<ji> runs over all pairs of 
nearest neighbor sites and the sum Y]^ ^ runs over all distinct pair of sites of the lattice; r»j is the distance (in crystal 
units) between sites i and j; S represents the ratio between the exchange Jo and dipolar Jd coupling parameters, 
where the energy is mesured in units of Jd, which is assumed always antiferromagnetic (Jd > 0). Hence, S > means 
ferromagnetic exchange copuling. 

There are few numerical results concerning the equilibrium statistical mechanics, i.e. the finite temperature phase 
diagram of this model. Maclsaac and coauthorsa have shown that the ground state of Hamiltonian (Q) is the antiferro- 
magnetic state for i5 < 0.85. For <5 > 0.85 the antiferromagnetic state becomes unstable with respect to the formation 
of striped domain structures, that is, to state configurations with spins aligned along a particular axis forming a 
ferromagnetic strip of constant width h, so that spins in adjacent strips are anti-aligned, forming a super lattice in 
the direction perpendicular to the strips. They also showed that striped states of increasingly higher thickness h 
becomes more stable as S increases from S = 0.85. Moreover, they showed that the striped states are also more stable 
than the ferromagnetic one for arbitrary large values of S, suggesting such a fliiase to be the ground state of the model 
for S > 0.85. Monte Carlo calculations on finite lattices at low temperatureo'0 gave further support to this proposal, 
at least for intermediate values of <5. Furthermore, such simulations have shown that striped phases of increasingly 
higher values of h may become thermodynamicallju. stable at finite temperatures for intermediate values of 6. This 
results are in agreement with other analytic onesutl. For small values of 5 the system presents an antiferromagnetic 
phase at low temperatures. At high temperatures, of course, the system becomes paramagnetic. 
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The dynamics of the model in the striped region is characterized by the formation and growth of magnetic domains, 
dominated by the competition between the exchange and the dipolar interactions. Monte Carlo studies of the dynamics 
at low temperatures have shown the existence of two different dynamical regimes, according to the value of 6. First, for 
6 > S c ~ 2.7 the magnetization relaxes exponentiallyu, with a relaxation time that depends both on the temperature 
and 5. For 5 < S c the magnetization presents a power law decay, with an exponent independent of S. Second, strong 
hysteresis effects appearau for <5 > (L, which are almost absent for 5 < S c . Finally, different types of aging behaviors 
have been observed in both regimesu. 

Aging effects, that is, history-dependence in the time evolution of correlations and response functions after the 
system has been quenched into some non-equilibrium state, appear, in a variety of ordered and disordered systems 
which are essentially out of equilibrium on experimental time scalesu. Aging can be observed in real systems through 
different experiments. A typical example is the zero-field-coolingii3 experiment, in which the sample is cooled in zero 
field to a sub-critical temperature at time to- After a waiting time t w a small constant magnetic field is applied and 
subsequently the time evolution of the magnetization is recorded. It is then observed that the longer the waiting time 
t w the slower the relaxation. 

Although aging can be detected through several time-dependent quantities, a straightforward way to establish it in 
a numerical simulation is to calculate the spin autocorrelation function 

C{t + t w ,t w ) = —J2(a i (t + t w )a i {t w )) (2) 

i 

where < ■ • ■ > means an average over different realizations of the thermal noise and t w is the waiting time, measured 
from some quenching time to = 0. 

A second quantity of interest is the conjugated response function to an external magnetic field hi(t): 

o,,,, ■ x _ 1 ^ d{a t (t + t w )) 

R(t + t w ,t w )--^ . (3) 

In a variety of disordered systems and also in some domain growth processes in ordered ones C(t + t w ,t w ) and RU^A- 
t w ,t w ) are found to satisfy the generalized fluctuation- disspation relation proposed by Cugliandolo and Kurchanliil: 

, \ X(t + t w ,t w ) dC(t + t w ,t w ) 
K(t + t w ,t w ) = - — (4) 

At equilibrium C(t + t w , t w ) and R(t + t w , t w ) satisfy time translational invariance (TTI), the functions depend only 
on the times difference t and X(t + t w ,t w ) = I, that is, Eq.(|^) reduces to the usual fluctuation-disspation theorem 
(FDT) . Out of equilibrium such properties are not expected to hold depending on the observation time scales. The 
following scenario has been proposed in the context of spin glasseaHl: for small values of t (t/t w << I) the system is 
in quasi-equilibrium and equilibrium properties hold; in the aging regime t/t w >> 1 both TTI and FDT do not hold, 
i.e., C(t + t w ,t w ) depends explicitly on t and t w and X (t + t w , t w ) ^ 1. Moreover, for large values of t w , X(t + t w ,t w ) 
becomes a function of time only through C(t + t w , i—l: X (t + t w , t w ) = X(G(t + t w , t w )). This function X(C) has 
been interpreted in terms of an effective temperaturea. 

At high temperatures X equals one since the system always equilibrate at large times and the equilibrium properties 
hold. At low temperatures, where aging phenomena appear, the departure of X(C) from I characterizes the FDT 
violation. _ _ 

ThiS|-scenario has been verified in sever^ models of spin glassesE3, in the Lennard- Jones glasstj, in kinetic Ising 
modelailand in polymers in random mediaMI. It has also been verified in the domain growth dynamics of ferromagnetic 
systemsEZI in dimensions d — 2 and 3, where X has been found to be zero in the aging regime. 

Instead of analyzing the response R(t + t w , t w ) we look at the integrated response function (proportional to the mag- 
netic susceptibility), that is, in a zero-field cooling numerical experiment we observe the growth of the magnetization 
under a constant external field applied at t w : 

rt+t w 

M(t + t w ,t w )= R(t + t w ,s)h(s)ds. (5) 



Using Eq.([|) we can rewrite Eq.(|]) for long times as 



^M{t + t w ,t w ) = I X(C)dC. (6) 

C(t+t m ,t m ) 



h 
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If FDT is satisfied Eq.(^|) reduces to a linear relation 

-M(t + t w ,t w ) = l-C(t + t w ,t w ), (7) 

while a departure from this straight line in an M vs C parametric plot indicates a violation of FDT and gives 
information about X(C). 

In this work we present the results of Monte Carlo simulations in the two-dimensional Ising model defined by the 
Hamiltonian (|l|) on an N = 30 x 30 square lattice with free boundary conditions. We chose the heat-bath algorithm 
for the spin dynamics and time is measured in Monte Carlo steps per site. For each run the system is initialized in a 
random initial configuration corresponding to a quenching from infinite temperature to the temperature T at which 
the simulation is done. We compute C(t + t w , t w ) as a function of the observation time t, for different values of t w , 8 
and T. In the striped phase at low temperatures this function decays quickly from C(t w ,t w ) = 1 to a constant value 
that jpersists for t/t w « 1; at t > t w it decays more slowly towards zero with a scaling law that depends on the 
ration h(t)/h(t w ). The scaling function h(t) appears to be linear for S > S c and logarithmic for S < S c . This behavior 
suggests a different domain growth dynamics in every one of the dynamical regimesu. Hence, it is of interest to check 
whether such different dynamics are reflected in the FDT violation factor or not. 

At time t w we take a copy of the_system, to which a random magnetic field h(i) = h ej is applied, in order to avoid 
favouring one the different phasesO; are taken from a bimodal distribution (e, = ±1) and the stresaht h of the 
field is taken small (h = 0.01) to ensure linear response. We then compute the staggered magnetizationliJ 

M (t + t w ,t w ) = —^2 (°i(t + t w )ti) (8) 

i 

whose conjugate field is h and where the overline means an average over the random variables e^. We then obtain 
parametric plots of T M(t + t w ,t w )/h vs. C(t + t w ,t w ) for several values of the temperature, the waiting time t w = 2™ 
(n = 5, 7, 9, 11) and for 6 — 2 < S c and 5 = 4 > S c . We first made some checks at high temperatures where the system 
equilibrates quickly, verifying that both TTI and FDT are satisfied. 

In Fig. (Jl|) a parametric plot of the integrated response vs. autocorrelation is shown for S = 2, T = 0.5 and waiting 
times t w — 2™ (n — 5, 7, 9, 11) from top to bottom. We made an average over 400 realizations of the random field. The 
straight line corresponds to the FDT relation ((7J) with constant slope —1. For fixed t w and observation times t «t w 
FDT holds, the system is in the stationary regime. For times t ~ t w the curves begin to depart from the FDT line 
signalling a crossover region where the system begins to fall out of equilibrium. Finally when t >> t w the system is 
out of equilibrium with the correlations decaying to zero as t — > 00. In this last regime the integrated response keeps 
growing for the small t w but, as t w grows it tends to stabilize in a constant value. Furthermore, this value decreases 
as t w grows. The older the system the smaller the memory of the past history. This behavior is similar, e.g. to 
what happens in the coarsening dynamics of ferromagnetalll. It is not at all obvious that this would be the case. It 
is important to note thatpthe system is in a region (<5 = 2) where the ground state is the striped phase with stripe 
width h which, grows witha 5. For this value of S the scaling of the autocorrelations in the out of equilibrium regime 
is logarithmiccl: C(t + t w , t w ) oc log(t)/log(r(t w )). This slow decay is typical of activated dynamics in systems with 
a broad distribuition of relaxation times. This may be consequence of the degeneracy of the striped ground state, a 
problem that deserves further study. So naively one would expect rather strong memory effects as a consequence of 
the slow logarithmic decay in the correlations, at variance with what is observed in the simulations. 
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FIG. 1. Integrated response versus autocorrelations for T = 0.5 and 8 = 2. The different curves correspond, from top to 
bottom, to waiting times t w = 2 5 , 2 7 , 2 9 , 2 11 . The straight line corresponds to the FDT relation TM/h = 1 — C. 

The flatness of the integrated response for long times implies an FDT violation factor X(t + t w ,t w ) = for t w — > oo. 
If we interpret TJX as the "effective temperature" for the system in this time regime, this implies an infinite effective 
temperaturet3. 

The FDT plot for 5 = 4 is presented in Fig.(||). This is qualitatevely similar to the plot for 6 = 2. The main 
difference is that the integrated response flattens to a value which is roughly half of that corresponding to 6 = 2. In 
this case the ferromagnetic term of the Hamiltonian is clearly dominant. This is reflected, e.g. in the scaling form of 
the autocorrelations in the aging regimeo, i.e. C(t + t w ,t w ) oc t/r(t w ). We must note, however, that the stable phase 
still corresponds to the striped oneH but with increasing value of the width of the stripes as 6 increases. For a fixed 
t w in the aging regime, the striped domains are wider with 5 = 4 than with 6 = 2. Consequently, the domain walls 
have smaller total lenght in the latter case. This implies a smaller contribution for the staggered magnetization while 
the bulk contributions are roughly the same. So it is reasonable to expect a smaller response in this case. In other 
words, as the systems becomes more "ferromagnetic" the long term memory becomes weaker. 
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FIG. 2. Integrated response versus autocorrelations for T 




The different curves correspond, from top to 



bottom, to waiting times t w = 2 a , 2 , 2 , 2 . The straight line corresponds to the FDT relation TM/h = 1 — C. 



In Fig.(||) we show a parametric plot for different values of the temperature T = 1, 0.7 and 0.5 (from top to 
bottom) and a fixed t w = 2 7 for 6 = 2. Equilibrium dynamics is restored at T = 1. It would be interesting to know if 
the change of dynamical regime at this temperature coincides with a thermodynamic phase transition or is a purely 
dynamic effect. 






FIG. 3. Integrated response versus autocorrelations for t w = 2 7 and 6 = 2. The different curves correspond, from top to 
bottom, to temperatures T = 1, 0.7, 0.5. The straight line corresponds to the FDT relation TM/h = 1 — C. 
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We have studied the FDT violation in the coarsening process of an Ising model with dipolar long range interactions. 
Going through 8 C ~ 2.7 the aging dynamics of the autocorrelations presents a crossover from a logarithmic decay for 
5 < S c to an algebraic decay for 5 > S c , probably related to some dynamical phase transition related with the change 
of the strips width. We asked weather this difference would manifest itself in the responses and FDT violation factor 
X . It turns out that this is not the case. For long waiting times in the aging regime X — > in both cases, signalling 
that the long term memory is weak in both regimes. 
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